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Abstract 

J> Because of the fine-tuning problem in classical Majoron models in recent years 

O several new models were invented. It is pointed out that double beta decays with 

Majoron emission depend on new matrix elements, which have not been considered 

in the literature. A calculation of these matrix elements and phase space integrals 

is presented. We find that for new Majoron models extremly small decay rates are 
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^ In many theories of physics beyond the standard model neutrinoless double 

beta decays can occur with the emission of new bosons, so-called Majorons 
[1-4]: 

2n ^ 2p + 2e- + </) (1) 

2n ^ 2p + 2e- + 2</). (2) 

Since classical Majoron models [1,5] require severe fine-tuning in order to 
preserve existing bounds on neutrino masses and at the same time get an 
observable rate for Majoron emitting double beta decays in recent years several 
new Majoron models have been constructed [6-8], where the terminus Majoron 
means in a more general sense light or massless bosons with couplings to 
neutrinos. The main novel features of the "New Majorons" are that they 
can carry units of leptonic charge, that there can be Majorons which are 
no Goldstone bosons [6] and that decays with the emission of two Majorons 
[4,7] can occur. The latter can be scalar-mediated or fermion-mediated. In 
vector Majoron models the Majoron becomes the longitudinal component of 
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Table 1 

Different Majoron models according to Bamert/Burgess/Mohapatra^. The case IIF 
corresponds to the model of Carone^°. 

a massive gauge boson [8] emitted in double beta processes. For simplicity we 
will call it Majoron, too. 

In tab. 1 the nine Majoron models we considered are summarized. [7,8] It is 
divided in the sections I for lepton number breaking and II for lepton number 
conserving models. The table shows also whether the corresponding double 
beta decay is accompanied by the emission of one or two Majorons. 
The next three entries list the main features of the models: The third column 
lists whether the Majoron is a Goldstone boson or not (or a gauge boson in 
case of vector Majorons IIF). In column four the leptonic charge L is given. 
In column five the "spectral index" n of the sum energy of the emitted elec- 
trons is listed, which is defined from the phase space of the emitted particles, 
G ~ {Qf5f5 — r)", where Qpp is the energy release of the decay and T the sum 
energy of the two electrons. The different shapes can be used to discriminate 
the different decay modes from each other and the double beta decay with 
emission of two neutrinos. In the last column we listed the nuclear matrix 
elements which will be defined in more detail later. Nuclear matrix elements 
are necessary to convert half lives (or limits thereof) into values for the effect- 
ive Majoron-neutrino coupling constant, using the approximate (see below) 
relations: [4,9] 

[Tii2]-^ = \<ga>r ■\M„\^ -Gbb^ (3) 

with m = 2 for /3/3(/)-decays or m = 4 for /3/3(/)(/)-decays. The index a in eq. 
(3) indicates that effective coupling constants (/„, nuclear matrix elements M„ 
and phase spaces GbBc differ for different models. 

As shown in tab. 1, several Majoron models with different theoretical mo- 
tivation can lead to signals in double beta decays which are experimentally 
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indistinguishable. The interpretation of experimental half life limits in terms 
of the "effective Majoron-neutrino coupling constant" is therefore model de- 
pendent. Subsequently we give a brief summary of the theoretical background 
on which our conclusions on the different Majoron models are based. 
Single Majoron emitting double beta decays (OufSfScj)) can be roughly divided 
into two classes, n = 1 (case IB, IC and IIB) and n = 3 (IIC and IIF) decays. 
As has been noted in [7] as long as OufSfS decay has not been observed, the three 
n = 1 decays are indistinguishable from each other. We will call these Majorons 
"ordinary", since they contain the subgroup IC, which leads to the classical 
Majoron models. [1,2,10] For all ordinary Majorons the effective Majoron- 
neutrino interaction Lagrangian, leading to OufSfScj) decay is [2,6] 

^?;f • = -l^a.^PL + h,PR)u,ct>* + h.c. (4) 

Here, Prjl = l/2(f i 75). Using eq. (4) the amplitude corresponding to the 
Feynman graph is, in the notation of [6] 

Vei^ Vej are elements of the neutrino mixing matrix, nii and nij denote neutrino 
mass eigenvalues and wpjaj are nuclear matrix elements containing double 
Fermi and Gamow-Teller operators. To arrive at the factorized decay rate 
eq. (3), the usual assumption <^ q ^ pp ^ O{100MeV)^ where pp is 
the typical Fermi momentum of nucleons, is made. By this assumption the 
term proportional to can be dropped and the effective coupling constant 
is defined as: 

{gf-M. = J2v^^V^,k,. (6) 

hi 

In this approximation matrix elements for ordinary Majoron decays coincide 
with the leading terms Mqt and Mp of the well-known mass mechanism of 
OufSfS decay. 

Burgess and Cline advocated the so-called charged Majoron model IIC. [6] In 
this model the effective interaction Lagrangian is 

L^f = -^uriALPL + AnPn)ud,^ + h.c. (7) 

Note, that in the charged Majoron model the two additional powers of n in the 
phase space integrals originate from the derivative coupling of the Majoron in 
L^'^' . As shown in [6], for charged Majorons the contribution from the leading 
order matrix elements to the decay rate vanishes identically, so that one has 
to go to the next higher order in the non-relativistic impulse approximation 
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of hadronic currents. The amphtude for OufSfScj) decay is then given by 

^■-'■(o.m) = 8v^E J 1^^-,— + ..,(8, 

which leads to an effective couphng constant {g)^'^' as in the ordinary Ma- 
joron case, but with hi^ given by 6 = j{Ai,m* + m*AR), with the neutrino 
mass matrix m, generator matrices Aj^jji and the decay constant /. The had- 
ronic term Wq is similar but not identical to the recoil matrix element of OufSfS 
decay induced by right-handed currents. This difference has turned out to be 
important. In the notation of [6] 

= I I,, 2^^^ , ■ A F\e-''%l{Cnal, - Oman) + (/^(D„ - D„)]|/),(9) 

= I i^/?^^^, ■ A F\e-''^-'[B^ X a- + D„ X (10) 

in which the summation over Y,nm '^n^m suppressed. Here, C„ and D„ are 
nuclear recoil terms [9] 

= (P. + P;) • a.J{2M.,) - {E^ - E'JiP^ - P'J ■ aj{2ml) (If) 

D„ = [(P„ + P;) + z/i^(P„ - P;) X a-;]/(2M„) (12) 

P„ (En) and P^ {E'n) are momenta (energies) of initial and final state nucleons, 
is the pion and M„ the nucleon mass and fip originates from the weak 
magnetism. 

The terms of are neglected compared to Wq due to the estimation (P„ + 
P;) < (P„ - K),{En - < 0{Qpp). [9] Following [II] we will also keep 
only the central part of the recoil term D. Although both are approximations, 
which needs to be checked numerically, we do not expect it to affect any of 
our conclusions. 

Finally for vector Majoron models (case IIF) [8] 

Llf = -^^Y{c,,Pl + d,,Pn)vX^ + h.c. (13) 
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where is the emitted massive gauge boson. The effective coupling constant 
can be defined as in the ordinary Majoron model, with the replacement hi^ = 
-^{cijmj — rjiidij), where M is the gauge boson mass. As discussed in [8], the 
vector Majoron amplitude approaches the charged Majoron one in the limit of 
vanishing gauge boson masses, which we assume in the phase space integration. 
They depend on the same nuclear matrix elements than the charged Majoron 
discussed above. We will therefore not repeat the definitions here. 
Double Majoron emitting decays (Oi//3/3</)</)), mediated by fermions, can have 
either spectral index n = 7 or n = 3, depending on whether the Majoron 
couples derivatively suppressed or not. [7] 



4 



In addition, in principle OufSfScfxj) decays could also be mediated by exotic 
scalars. The amplitude of scalar-mediated decays, however, is expected to be 
very much suppressed, since the scalars must have masses larger than about 
50 GeV due to the LEP-measurements. [7] We will therefore concentrate on 
the fermion-mediated decays. 

The Yukawa coupling of the Majoron to the neutrinos for the n = 3 decays 
(cases ID, IE, IID) is given as: 

Lit = -HAaPL + B,aPR)Na<P + h.C. (14) 

where Aia and Bia represent arbitrary Yukawa-coupling matrices and Na are 
sterile neutrinos. The corresponding amplitude for OufSfScfxj) decay is 



[q^ — m? + ie){q^ — m| + ie){q^ — ml + ie 



Although for A^'^'{0iyf3f3(f)(f)) the same combination of nuclear operators ap- 
pears (wp — wgt), note the additional (g^ — m^)~^ compared to A'^'^'{0iyf3f3(f)). 
jVija in (15) is given by 

J^zja = -q^{A,aBjam^^ + Aj^i^^mj,^ + B.aBjarriNa) + A.aAjam^m^mNai^^) 

In order to separate the particle physics parameters from the nuclear structure 
calculation, it is most convenient to neglect the last term in eq. (16). This can 
be justified by considering that the mass eigenvalues rrii,^^ <^ pp so that the 
last term in eq. (16) for not too large mjVa is suppressed compared to the 
first three by at least m^^^/pp ~ (9(10~^~^). Then, the q^ is absorbed into 
the neutrino potentials and we redefine Afija to obtain the effective coupling 
constant as 

(g) = (— Yl Ve^Vej[A,aBjam^,^ + AjaB^aTTl^,. + B.aBjamNa])^ (17) 

nip ^ 

Note, that we have arbitrarily absorbed a factor of into the definition of 
((/) here to get for the effective coupling a dimensionless quantity. 
Eor the n = 7 OufSfScfxj) decays, the effective Lagrangian is (IIE/fermion medi- 
ated): 

L^.. = -tlAYiX.aPL + Y,aPR)Nad,cf> + h.C. (18) 

Again, Na denotes a sterile neutrino and the derivative coupling of (j) accounts 
for the additional powers of n in the phase space integrals. The amplitude for 
n = 7 decays is the same as for the n = 3 case, discussed above, with the 
replacement: J\f,ja = X.aY.aniNa- 
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Note, that Xia and Yia have the dimension of an inverse mass. Therefore, also 
((/) has a dimension of an inverse mass. To define a dimensionless couphng 
constant in this case one would have to specify the symmetry breaking scale, 
which is however undetermined by the model. 

For the Majoron models considered in this work there are five nuclear matrix 
elements to be calculated. Within the closure approximation they are defined 
as: 

Mf = (4) < Nf\\h.r,ass{l^,r)T + T+\\N, > (19) 

9a 

Mgt =< Nf\\h^ass{l^,r)T+T+ana^\\N, > (20) 

McR = (— )^ < Nj\\hn{fi,^T+T+a^a^\\N, > (21) 
9a 3 

Mf.^ = (4) < Nf\K2{i,,r)T+T+\\N, > (22) 
9a 

=< Nf\\h^2{ij,r)T+T+a.,a.r,\\N, > (23) 
where ha denote the neutrino potentials 



R r d^q e'"~ 
27r^ J to to -\- jjL 



hma.s.sij^,r) = — / _ _ , (24) 



Here /i = {En — Ej) denotes the average excitation energy of the intermediate 
nuclear states, to = + is the energy of the neutrino and since we 
assume all neutrinos to be light, the indices on neutrino masses have been 
dropped. Note, that in order to define matrix elements dimensionless we follow 
the convention of [9]. That is hmassi^^) and h^2[f) are arbitrarily multiplied 
by the nuclear radius R = VqA' with Tq = 1.2 fm, while hR^r) includes the 
nucleon mass. Compensating factors appear in the pref actors of the phase 
space integrals. 

We have carried out a numerical calculation of these matrix elements within 
the pn-QRPA model of [12,13]. To estimate the uncertainties of the nuclear 
structure matrix elements the parameter dependence of the numerical results 
has been investigated. Since the matrix elements Mgt and Mf have been 
studied before, [12] we will concentrate on Mcr, Mqj^2 and Mf,^2. Mgt and 
Mf can be calculated with an accuracy of about a factor of 2[12]. 
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The matrix element Mcr shows a very similar behaviour as Mqt- This is 
in agreement with the expectation, since only the central part of the recoil 
terms is taken into account, so that apart from the different neutrino potential 
Mcr has the same structure as Mqt- Neither variations of the strength of 
the particle-particle force gpp nor a change in the intermediate state energies 
significantly affects the numerical value of Mcr- We therefore conclude that 
Mcr should be accurate up to a factor of 2, as is expected for Mqt- 
Unfortunately, in the case of the matrix elements Mqj^2 and Mp^2 the situ- 
ation is very different. Both, variations of g^p or /i, can change the numerical 
results drastically (fig. f). In fact, it is found that Mqx^2 displays a very sim- 
ilar dependence on g^p as has been reported in pn-QRPA studies of 2i//3/3 
decay matrix elements. [12] Especially important is that in the region of the 
most probable value of gpp Mqj^2 crosses zero. 

Also for variations of the assumed average intermediate state energy a rather 
strong dependence of the results on the adopted value of [i has been found. As 
a consequence of this unpleasant strong dependence, for an accurate prediction 
of Mqj^2 and Mp^2 it seems necessary to go beyond the closure approxima- 
tion. 

The basic reason for the unusual sensitivity of Mqj^2 and Mp^2 on [i can be 
traced back to a certain difference in the neutrino potential of these matrix 
elements compared to M^t/f, hmassifJ-^'f^ ~ while /i^2(/i,r) ~ cu"'*. Con- 
tributions from very low momenta are therefore much preferred in /i^2(/i,r) 
compared to /imass(/W, r). (Note, that this leads also to a much smaller value 
for a typical uj than the naive expectation of ~ ~ (9(50 — 100) MeV!). 
With typical uj of only (9(few) MeV the strong dependence of /i^2(/i,r) on /i 
becomes obvious. 

Results of the calculation for various experimentally interesting isotopes are 
summarized in table 2. Note that the matrix elements are valid for the limit 
of small intermediate particle masses, up to the order of 10 MeV. If any of the 
virtual particles in the Feynman graphs can have masses larger than 10 MeV, 
the matrix elements are no longer constant and the values in table 2 should 
only be taken as upper limits for the analysis of data. 

In comparison to the nuclear matrix elements phase space integrals can be 
calculated very accurately, so uncertainties of this calculation will not be dis- 
cussed. We define the phase space integral as 

GbBc, = aa- {Qf3f3 - ei - £2)" n Pk<^kf{<^k)dek (27) 

k 

where the prefactor a„ depends on the Majoron mode under consideration. 
A summary of the definitions is given in table 3. Qpp is the maximum decay 
energy, Ck and pk are the energies and momenta of the outgoing electrons and 
f{€-k) is the Fermi function calculated according to the description of. [9] Note 
the large difference in the phase space values of the old (n = I) and new 
Majoron models. 



7 



nucleus 


Mp - Mr<T 






Co 7R 

yjG 1 D 


4.00 


U.iD 


^ in-3±i 




4.Uo 


n 1 /I 

U.i4 


^ in-3±i 


i\ j„ 1 nn 

iViO iUU 


/I 8K 
4.oD 


U.iD 


^ in-3±i 




Q 9Q 


n 1 n 

U.iU 


^ in-3±i 


To 1 98 


zL zLQ 




~ in-3±i 

~ lU 


Te 130 


3.90 


0.12 


~ 10-3±1 


Xe 136 


1.82 


0.05 


~ 10-3±1 


Nd 150 


5.29 


0.15 


~ 10-3±1 



Table 2 

Dimensionless nuclear matrix elements of Majoron emitting modes calculated in 
this work 
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Table 3 



Values of phase space integrals calculated in this work 

Having calculated nuclear matrix elements and phase space integrals, it is 

straightforward to derive limits on the effective Majoron-neutrino coupling 

constants for the various Majoron models from experiment. 

Although experimental half life limits are comparable for all decay modes, 

as observed recently for ''^Ge decay [14,15], restrictive limits on the coupling 

constants of ordinary Majoron models contrast with limits on any of the new 

Majoron models, which will be weaker by (3-4) orders of magnitude. 

The surprisingly weak limits which one obtains for the neutrino-Majoron 

coupling constant due to small matrix elements and phase spaces for all of 

the new Majoron models, require further explanation. (Note that the foUow- 
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ing discussion is independent of the isotope under consideration.) Consider, 
for example, ordinary and charged Majoron OufSfScj) decays. Limits on the 
effective couphng constant for single Majoron emitting decays will scale as 
{g) ~ A~^{Tif2 • Gbb)~^ ■ Thus, the relative sensitivity of a double beta de- 
cay experiment on ordinary and charged Majoron decays can be expressed as 

1 /2 

Inserting the definitions of the corresponding amplitudes, it is clear that even 
if the half life limit derived for the charged Majoron decay equals that of the 
ordinary Majoron mode, limits on the charged Majoron-neutrino coupling 
constant will be weaker by Mn/iQpp — T) ~ 1000 ! (Note, that this crude 
estimation is to first approximation independent of nuclear structure proper- 
ties.) 

A similar analysis can be easily done for double Majoron emitting decays. 
Again, very crudely, a reduced sensitivity of (487r^) • Pf l{Q(5(5 — T) ~ (few) 
xIO"* for n = 3 double Majoron decay, compared to ordinary Majoron decays, 
is expected. Here, the factor (487r^) is due to the phase space integration 
over the additional emitted particle, while the latter factor comes from the 
additional propagator. 

One might think that since our definition of the effective coupling constant for 
the n = 3 OufSfScfxj) decays includes a factor mN^/nie^ where mjVa is the sterile 
neutrino mass, one could get {g) easily as large as wanted, since the mass of 
the sterile neutrino is not bounded experimentally. However, matrix elements 
will fall off M ~ as soon as rriN^ is larger than the typical momenta. 

While for the matrix elements Mq^jp for ordinary Majoron decays such a 
reduction occurs starting from masses of exchanged virtual particles in the 
region of 100 — 1000 MeV, for Mqx^2 ip^2 the suppression will be important 
already for much smaller masses (see the E'n-dependance fig. I). 
Since the sensitivity of double beta decay experiments to the new Majoron 
models is so weak, it might be interesting to compare expected half lives for the 
different models for different ((/), ((/) lO""* as a typical sensitivity in coupling 
constant for ordinary Majoron models and ((/) = I as an upper possible limit 
allowed by perturbation theory, with current experimental limits of (9(10^^) 
years (see tab. 4). From this consideration it is very unlikely that any of the 
new Majoron models can produce an observable rate in planned or ongoing 
double beta decay experiments. Only the charged and the vector Majoron 
model [6,8] could produce an observable effect if Xlij KiKj is not smaller than 
O.I and the real coupling constant of order 0{1). 
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